Four symmetry classes of plane partitions under one roof 
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In previous paper, the author applied the permanent-determinant method of Kasteleyn and its non-bipartite 
generalization, the Hafnian-Pfaffian method, to obtain a determinant or a Pfaffian that enumerates each of the ten 
symmetry classes of plane partitions. After a cosmetic generalization of the Kasteleyn method, we identify the 
matrices in the four determinantal cases (plain plane partitions, cyclically symmetric plane partitions, transpose- 
complement plane partitions, and the intersection of the last two types) in the representation theory of sl(2,C). 
The result is a unified proof of the four enumerations. 



Stanley JlO[ | and Robbins recognized that there are ten sym- 
metry classes of plane partitions in a box, and, with the aid 
of computer experiments and other work, that the number of 
plane partitions in each symmetry class is given by a prod- 
uct formula. The program of proving each of these formulas 
has recently been completed [j], Bl but although the for- 
mulas are very similar, there is no known unified treatment. 
Moreover, some of the enumerations presently require diffi- 
cult, ad-hoc calculations with generating functions or matri- 
ces. In fact, the enumeration of one of the classes (totally sym- 
metric plane partitions or TSPP's) admits a natural ^-analogue 
which is still open. 

In a previous paper, the author [||| discussed the permanent- 
determinant method of Kasteleyn and its generalization, the 
Hafnian-Pfaffian method, as the first step towards a possible 
unified treatment. (The connection between this method and 
plane partitions was discovered jointly with James Propp.) 
In this paper, we carry out this strategy for the four sym- 
metry class that are given by determinants and not merely 
Pfaffians: Unrestricted plane partitions (PP's), cyclically sym- 
metric plane partitions (CSPP's), transpose-complement plane 
partitions (TCPP's), and cyclically symmetric, transpose- 
complement plane partitions (CSTCPP's). The new idea is to 
find the matrices given by the permanent-determinant method 
in the representation theory of the Lie algebra sl(2,C). The 
main result of the present paper is the following theorem: 

Theorem 0.1. A Kasteleyn-flat, weighted adjacency matrix 
for the graph Z(a,b,c), whose matchings are bijective with 
plane partitions in an a x b x c box, arises as a(X)\_ l 
in a representation a of sl(2,C) which is a tensor prod- 
uct of three irreducible representations, with a basis formed 
from weight bases of irreducible representations. Conse- 
quently, (detff(X)|_j)//n is the number of such plane parti- 
tions, where m is the value of any term in the determinant. 

Interestingly, the representation theory of the quantum 
group C7 ? (sl(2,C)) similarly yields the ^-enumeration of PP's, 
but the author could not obtain the known ^-enumeration of 
CSPP's by this method. 

Proctor's minuscule method is a competing method that, 
in various forms, yields the enumeration of five symmetry 
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classes of plane partitions [Q, |8, 11 1. Proctor's method partly 
explains Stembridge's q = — 1 phenomenon jl2|], and it is ap- 
plicable to any symmetry class not involving cyclic symmetry. 
Likewise, the method given here is applicable to any symme- 
try class which does not involve transposition or complemen- 
tation, but which may involve their product. 

The author believes that one of these two methods, or per- 
haps some combination of the two, should lead to a complete 
unified enumeration. Both methods involve Cartan-Weyl rep- 
resentation theory. Moreover, many of the enumerations can 
be analyzed with the Gessel-Viennot non-intersecting lattice- 
path method [||], which in turn can be recognized as a con- 
densed version of the permanent-determinant method [|]| . For 
example, the Gessel-Viennot method applied to PP's yield 
Carlitz matrices. The determinants of these matrices are well- 
known; interestingly, Proctor Q| rederived them, and there- 
fore enumerated PP's, by an argument which is parallel to the 
one presented here. 

The author would like to thank James Propp for suggest- 
ing numerous corrections and for his continued interest in this 
work. 



1. PRELIMINARIES 

1.1. Plane partitions and their symmetries 

A plane partition in an a x b x c box is a collection of unit 
cubes in the rectangular solid [0, a] x [0, b] x [0, c] which is sta- 
ble under gravitational attraction towards the origin. Equiva- 
lently, it is an a x b matrix of integers with entries between 
and c; the stack of cubes is then a Ferrer's diagram (or 
bar graph) of the matrix. Figure [l] shows a picture of a plane 
partition together with three back walls of the box. We let 
N(a,b, c) be the total number of plane partitions in an a x b x c 
box. 

There are three natural symmetry operations on plane par- 
titions: 

1. Transposition, T, is a reflection through a diagonal plane of 

the box of a plane partition. A transposition-symmetric 
plane partition is called, briefly, a symmetric plane par- 
tition. 

2. Rotation, p, is a rotation by 120° about the long diagonal 

of the box of a plane partition. Symmetry under rotation 
is also called cyclic symmetry. 
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Figure 1 : A plane partition in a box and its lozenge tiling. 



3. Complementation, K, takes the point (x,y,z) to the point 
(a — x, b — y, c — z) for anaxbxc box. The operation K 
is the same for individual cubes, but for a plane partition 
P, k(P) is the set of all cubes C such that jc(C) ^ P. 

By convention, the operations T and p are only defined 
when the hexagon or box possesses them as symmetries, while 
the operation K always exists. The group of all symmetries has 
ten conjugacy classes of subgroups, and for a subgroup G, let 
N G (a,b,c) be the number of plane partitions with that symme- 
try. We describe G by generators as G = (g l ,g 2 ,...) and we 
often omit the angle brackets, so that N p r (a,a,a) is the num- 
ber of totally symmetric plane partitions. (By convention, the 
phrase "total symmetry" does not include the symmetry oper- 
ation of complementation.) 

As mentioned in the introduction, we will consider the four 
symmetry classes whose symmetry groups are subgroups of 
(p , JCT). In other words, CSPP's are invariant under p, TCPP's 
are invariant under kt, and CSTCPP's are invariant under 
both. 

In addition to ordinary enumeration, we will also con- 
sider the weighted enumeration of plane partitions, where 
the weight of a plane partition with n cubes is q" . Specifi- 
cally, let N G (a,b,c)q be the total weight of all plane partitions 
in the given symmetry class. The polynomials N(a,b,c)q, 
N x (a,a,b)q, and N p (a,a,a)q have nice product formulas, but 
we will only give a proof of the formula for N(a,b,c) q . 
(There is an alternate weighting whereby the weight of a 
plane partition is q" if it has n orbits; the corresponding poly- 
nomial N' t (a,a,b)q is known to have a product formula and 
N' p T (a,a,a) is conjectured to do so, but we will not consider 
these interesting cases.) 



1.2. Product Formulas 

We define some products related to Stanley's formulas Jic|]. 
If a { , a 2 , ■ ■ ■ , a k are positive integers and n is a non-negative 
integer, define the box product C(a l ,a 2 , ...,a k ;n) by: 



C(a v a 2 ,...,a k ;n) = ]~[ f~[ Yl max(n- J^x ( , 1) 
x l =0x 2 =O x k =0 i'=l 

A typical box product is the product of the numbers shown in 



Figure^ Define the simplex product T(k,n) by: 

k 

T(k,n) = C(°°,°<v . . ,°°;n). 
Thus, T(l,n) = n\, while 

T(2,n) =n!(n-l)!(n-2)!...2!l!. 
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Figure 2: Factors of the box product C(6,4;7). 

The simplex products can also be defined inductively by the 
following rules. 

T(0,n)=n 
T(k,0) = 1 

T(k,n)=T(k-l,n)T(k,n-l) (1) 

for positive integers k and n, with T(k,n) — 1 for all other 
cases. 

It follows immediately from the definition that 



C(flj,a 2 , • • • ,a,;n) 



C(°°,a 2 ,a v ...,a k ;n) 
C(°°,a 2 ,a^, . . . ,a k ;n — Oj) 



Applying this formula to each variable in turn yields a multi- 
plicative inclusion-exclusion formula: 



C(aj,a 2 • • • ,a k ;n) = 
T(k,n)T(k,n — a 
T(k,n^a l — a 2 )T(k,n — a l — « 3 



T(k,n)T(k,n — a l ) l T(k,n — a 2 ) 1 . . .T(k,n — a k ) 1 



In particular, 

C(a,b,c;n) = 
T(3,n)T(3,n — a — b)T(3,n — a — c)T(3,n — b — c) 
T(3,n — a)T(3,n — b)T(3,n — c)T(3,n — a — b — c) 

Similarly, define a box ^-product by: 



(2) 



C(a^,a 2 , . . . : a k ;ri/q 



n 



where (n) q = l+q + ... + q" , and a simplex g-product by 
T(k,n) q = C(o°,...,o ;n) q . 
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The formula for N(a,b,c) given by Stanley [ |lOp and origi- 
nally proved by MacMahon, is 



N(a,b,c) 



C(a,b,c;a + b + c — 1) 
C(a,b,c;a + b + c — 2) 



Combining this equation with equation ([|) and equation ([!]) 
we obtain the expression: an expression 



N(a,b,c) = 



T(2,d- \)T{2,a- l)T(2,b- l)T(2,c- 1) 



T(2,c-l)T(2,b-l)T(2,d-l) 
where, for abbreviation 



a = b + c 
c = a + b 



I ? = a + c 
d = a + b + c. 



This expression was suggested by Propp [||] in the form 

H(a + b + c)H{a)H{b)H{c) 



N(a,b,c) 



H(a + b)H(a + c)H(b + c) 



where 



ff(/i) = l!2!...(n-l)!=T(2,n-l) 



is the hyperfactorial function. 

For either formulation of MacMahon's enumeration, we 
can g-ify throughout. 



2. THE ENUMERATION 
2.1. The permanent-determinant method 

Following reference [0], a plane partition in an a x b x c 
box is equivalent to a tiling of a certain hexagon H(a,b,c) by 
unit lozenges, where a lozenge is a rhombus with a 60 degree 
angle, as shown in Figure |] The hexagonH(a,Z?,c) has angles 
of 120 degrees and edge lengths of a, b, c, a, b, and c, going 
clockwise around the perimeter. The hexagon has a unique 
tiling by unit equilateral triangles, and a unit lozenge tiling is 
therefore equivalent to a perfect matching of a bipartite graph 
Z(a,b,c), shown in Figure |[ whose vertices are the triangles 
and whose edges are given by adjacency. 

The number of perfect matchings of an arbitrary bipartite 
graph is given by the permanent of its bipartite adjacency 
matrix. Kasteleyn's permanent-determinant method dictates 
that, if the graph is planar, then there is a way to change 
the signs of the permanent to convert it to a determinant. 
(More accurately, Kasteleyn considered the more general non- 
bipartite case and produced a Pfaffian; the bipartite case was 
clarified by Percus ftfo.) Indeed, the permanent of the adja- 
cency matrix A(a,b,c) of Z(a,b,c) equals the determinant up 
to a global sign, and the global sign is in any case ambigu- 
ous because the rows and columns are unordered. One way to 
demonstrate the equality is to note that two matchings differ 
by a 3-cycle (carried by a hexagon of Z(a,b,c)) if their plane 
partitions differ by a cube. Since a 3-cycle is an even permuta- 
tion, it incurs no sign change between the corresponding terms 




Figure 3: The graph Z(2, 2, 3). 



of DetA(a,b,c). Since any two plane partitions are connected 
by the operation of adding or taking away an individual cube, 
it follows that all terms have the same sign in the determinant. 

More generally, suppose that Z(a,b,c) is weighted with 
non-zero weights, and that M is its weighted adjacency ma- 
trix. Then the non-vanishing terms in the determinant of M 
also correspond to matchings of Z(a,b,c), but they are no 
longer necessarily equal. If, for every hexagon of Z(a,b,c), 
the weights x l ,y l ,z l ,x 2 ,y 2 , andz 2 shown in Figure ^| satisfy 

then the non-vanishing terms of DetM are equal. In this case, 
we say that M and the weighting are Kasteleyn- flat. The num- 
ber of matchings is therefore 



(DetM)/m, 



(3) 



where m is the value of any single term. This principle can be 
extended to the ^-enumeration problem: Define the Kasteleyn 
curvature of a hexagon as 



(This statistic is called curvature by analogy with the curvature 
of a connection on a line bundle or the coboundary operation 
in homology; see reference |Q| for details.) If the Kasteleyn 
curvature is q everywhere and t is the weight of the matching 
of the empty plane partition, then 

N{a,b,c) q = (DetM)/m. 

In this case, the value of a term goes up by q if we add a cube 
to the corresponding plane partition. 

The symmetry operations t, p, and K act on hexagons 
H(a,b,c) respectively as rotation by 120°, rotation by 180°, 
and reflection about a diagonal. Therefore kt acts as reflec- 
tion about a bisector (a line which meets two opposite edges 
of the hexagon in the middle), and in its action on Z(a,a,2b), 
it fixes a row of edges which separates the graph into two 
isomorphic subgraphs. (Note that there are no K"T-invariant 
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Figure 4: The weights that determine Kasteleyn curvature. 



matchings of Z(a,a,2b — 1), or equivalently N KT (a,a,2b — 
1) =0.) The fixed edges must appear in any invariant sub- 
graph, so that TCPP's correspond to matchings of either sub- 
graph, denoted by Z KZ (a,a,2b). Similarly, (p,Kz) has three 
bisectors as lines of reflection, so that CSTCPP's correspond 
to matchings of any of the six subgraphs Z p KZ (2a,2a,2a) that 
remain. Since all the faces of these graphs are also hexagons, 
the same analysis of Kasteleyn's method apply fully, except 
that Kasteleyn curvature of q is unnatural from the point of 
view of plane partitions unless q — ± 1 . We will only consider 
q = 1 for these cases in this paper. 

A p-invariant matching of Z(a,a,a) is simply a matching 
of the quotient graph Z p (a,a,a) — Z(a,a,a)/p, and again the 
same analysis of Kasteleyn's method applies, except that the 
face in the center is now a 2-gon. This means that Z p (a,a,a) 
is not a simple graph, the matrix M is not uniquely determined 
by the weighting of Z p (a,a,a), and the matchings are not bi- 
jective with the non-zero terms in the determinant of M. How- 
ever, if the weight of each edge is a separate variable, then the 
monomials in the expansion of detM are bijective with the 
matchings of Z p (a, a, a). Taking this approach, if the 2-gon 
has sides with weight x { and x 2 , we define its Kasteleyn cur- 
vature as x 1 /x 2 and say that it is flat if x l = x 2 . As before, all 
terms in the determinant are equal if the weighting is flat at 
every face. More generally, the assignment of Kasteleyn cur- 
vature which corresponds to the ^-enumeration N p (a, a, a ) q is 
q 3 for every hexagon and q for the 2-gon. 



2.2. Representation theory of sl(2, C) 

In this section, we will find a Kasteleyn-flat matrix M for 
Z(a,b,c) in the representation theory of sl(2,C). We begin 
by reviewing the finite-dimensional representation theory of 

sl(2,C). 

The Lie algebra sl(2,C) of traceless 2x2 matrices has an 
important basis H, X, and Y given by the matrices: 



H = 
X = 
Y = 







1 





1 



extension of the Lie bracket on this basis: 

[H,X] = 2X 
[H,Y] = -2Y 
[X,Y]=H 

Let x and y be the standard basis vector on which the above 
matrices act. For each positive integer n, sl(2,C) has an ir- 
reducible representation V„ of dimension n + 1 which can be 
described as the action of sl(2,C) on homogeneous polyno- 
mials in x and y of degree n. Let a„ be the representation map 
from sl(2,C) to End(V„). To derive the action of sl(2,C) in 



the the monomial basis je",je" l y,x n 2 y 2 
nient to express X, Y, and H formally as: 

X =x^- 
dy 

y dx 

H = x-- - 
dx dy 

The corresponding matrices are then: 



,y", it is conve- 



a n (H) 



(n •• 
n-2 •• 
n-4 •• 



a„(X) 



\o 

/0 1 

2 

3 






-n) 



a n (Y) 





Vo o o o 




n-1 
n 








n 

0/ 



n 



\0 



0\ 





1 0/ 



The monomial basis is also known as the (dual) weight ba- 
sis, and we rename the basis vectors e n ,e n 2 , . . . ,e_ n , where 
the eigenvalue of e ; with respect to H is ;'. 

Theorem 2.1. (Clebsch-Gordan) Any finite-dimensional, ir- 
reducible representation of sl(2,C) is isomorphic to a direct 
sum of V„ 's. In particular, their tensor products decompose 
according to the equation: 



V n ®V k ^V n+k ®V n+k _ 2 ®V n+k _ A t 



\n-kY 



The entire Lie algebra structure can also be defined by linear 



We assume the first half of the theorem without proof. 
To understand the Clebsch-Gordan decomposition of a tensor 
product, we recall that if L is an element of an arbitrary Lie 
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algebra, its action on the tensor product of two representations 
a and j3 is given by: 

(o®/3)(L) = a(L)®/+/®j3(L), (4) 

where / is the identity matrix. 

Given the first half of the theorem, the following is an 
argument for the Clebsch-Gordan formula: Let S(H) and 
5(sl(2,C)) be the semi-ring of reducible representations of H 
with integer eigenvalues and the semi-ring of representations 
of sl(2,C), respectively. The elements of these semi -rings are 
representations and the algebraic structure is given by © and 
<g>. The semi-rings extend to Grothendieck rings R(H) and 
/?(sl(2,C)) by introducing subtraction. Since each V„ restricts 
to a reducible representation of H, 5(sl(2,C)) C S(H) and 
/?(sl(2,C)) C R(H). The ring R(H), as an abelian group, is 
generated by one-dimensional representations E n with eigen- 
value n for some integer n. By equation (Q), eigenvalues add 
under tensor product, and therefore there exists an isomor- 
phism ch : R(H) — > Z[f , 1 /t] given by the formula ch(£„) = t" . 
This isomorphism has a restriction ch : /?(sl(2,C)) — * Z[f , 1 /?] 
that we will call the character map. By the form of a„(H), 

ch(v n ) = t" + t n - 2 + ... + r". (5) 

Because the ch(V„)'s are linearly independent, ch is injec- 
tive. The Clebsch-Gordan formula then follows easily by ap- 
plying ch to both sides. Indeed, the assertion that the character 
map is injective on S(sl(2,C)) is a convenient restatement of 
the Clebsch-Gordan theorem. 

Since the matrices for a n (H), a„(X), and U n {Y) are nearly 
permutation matrices, they can be described graphically as in 
Figure |[ Each dot in the figure represents a basis vector which 
is an eigenvector of a n (H) with the given eigenvalue. The 
maps a„(X) and a„(Y) send each basis vector to another basis 
vector times the given scalar factor, except for the vector at 
the right (resp. left), which is in the kernel of a„(X) (resp. 
On{Y)). 



X 6 5 4 3 2 1 



H -6 -4 -2 2 4 6 



Figure 5: An irreducible representation of sl(2). 

As a warm-up for the main construction, consider the action 
of X, Y, and /Zona tensor product representation V = V n ® V, 
in the tensor product basis {ej(£>e-}. If a is the representation 
map, the actions of a(H) and tt(X) are represented diagram- 
matically by Figure ^. Each dot is a basis vector as before, 



4\ 3 2 1 
• • • ► • ► • 




• ► • ► • ► • ► • 




• ► • ► » N ► » N ► • 




a(X) V\_ x V\ 



Figure 6: The representation V 4 ® 



and diagonal strings span eigenspaces of tt(H) with the given 
eigenvectors, because e ; ® e ; - is an eigenvector of tt(H) with 
eigenvalue i + j. By equation (^), most columns of the ma- 
trix for a(X) have two terms, which are given by the arrows. 
(Technically speaking, one might say that the diagram is a 
directed, weighted graph and the matrix for (x(X) is the asym- 
metric weighted adjacency matrix of the graph.) Let VL be 
the eigenspace of u(H) with eigenvalue A. Then the matrix 
for a(X) can be divided into blocks a(X)\^ '-V\ x — > V|^ +2 
which are maps between adjacent eigenspaces. For example, 
in Figure |,a(X)|_|, whose domain and target are delineated, 
has the matrix: 

/4 1 0\ 

3 2 01 

2 3 
\0 1/ 

Consider the graph Z(a,b,c) as in Section [2.1[ let 

v = v e _ l ®v hl ®v a _ v 

and let a be the representation map. The tensor product ba- 
sis {e j (g) e j ® e k } for this representation can be depicted as a 
3 -dimensional rectangular box of points. Unfortunately, it is 
hard to show all of the points in a 2-dimensional picture. How- 
ever, a(H) is diagonal as before, and the basis vectors lying 
in each eigenspace constitute a 2-dimensional slice of the box, 
because each e i ® e ■ ® e k is an eigenvector with eigenvalue 
i + j + k. These slices are given as a sequence for V 4 £g> V 4 (g> V 5 
in Figure ^, together with the eigenvalue of a(H) for each 
slice. We will consider the slices V| i and V | , , which are 
given separately in Figure |[ together with points correspond- 
ing to their basis vectors. 

Most columns of the matrix for a(X) now have three 
terms, and the matrix divides into submatrices a(X)\* : VL — > 
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Figure 8: Basis vectors in V\_ l and V\ 



V\^ +2 . Consider the map a(X)\_ l , which is diagrammed in 
Figure ||. If we replace the arrows in this diagram by unori- 
ented edges, we see that a(X)\_ l is a weighted, bipartite ad- 
jacency matrix of the graph Z(a,b,c). Moreover, the weight- 
ing is Kasteleyn-fiat, because the weights of opposite edges 
of each face are equal. We have proved the following propo- 
sition: 

Proposition 2.2. If the representation V ±1 ® Vg { <S> V~_ l of 
sl(2,C) has representation map (X and is given with the basis 
which is the tensor product of the weight bases of the three 
factors, then the matrix for a(X)\_ l , the restriction of a(X) 
to the eigenspace of a(H) with eigenvalue — 1, is a Kasteleyn- 
fiat adjacency matrix for Z(a,b, c). 

Applying equation (||), 

N(a,b,c) = (det(a(X)\_ 1 ))/m, (6) 



where m is the value of any single term in the expansion of 
the determinant. Indeed, in this case all terms are not only 
equal, but are products of the same weights. The product of 
the weights of the horizontal edges in any matching is: 

a b 

Y\Yl(i + j-l) = C(a,b-a + b-l) 

i=U=l 

Therefore 

t = ±C(a,b;c- l)C(a,c;b- l)C{b,c;a- 1). (7) 

The map a(X)\_ l is a linear transformation between differ- 
ent vector spaces. For the purpose of computing determinants, 
it is easier to work with an endomorphism of a single vector 
space. Observe that the restriction a (F) | j :V\ l — >V|_] has all 
of the same properties as a(X)\_ l , and that the composition 
0C(F)| 1 0t(X)|_ 1 is an endomorphism of V | j . In particular, 

N{a,b,c) 2 =det(a{Y)a(X)\_ l )/m 2 . 

Since a(Y)a(X)\_ l comes from the action of sl(2,C) on 
V, it acts on each summand of a direct-sum decomposition 
of V separately. Therefore such a decomposition diagonalizes 
a(Y)a(X)\_ l and can be used to find its determinant. 

Lemma 2.3. The map 

is a map from a 1 -dimensional vector space to itself and its 
effect is multiplication by n 2 . 

Compare with Figure || or with the definition of a„ (X) and 
a„{Y). 
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In general if a and j3 are any two representation of sl(2, C), 

det((a + /3)(7)(a + J 3)(X)L 1 ) = 

det(a(Y)a(X)\_ l )det(l5(Y)l5(X)\_ l ). 

It follows that the assignment 

a ^ det(a(Y)a(X)\_ l ) 

defines an abelian-group homomorphism from the additive 
subgroup of the Grothendieck ring /?(sl(2,C)) spanned by 
V 2 j's to the multiplicative group Q* of non-zero rationals. 
This homomorphism factors through the character map ch to 
yield a map D such that 

D(ch(V))=det(a(F)a(X))|_ 1 ) 

for any odd-weight representation V. By equation (||), we can 
define D by 



D(t 2 "- 1 ) 



{n-lf 



for n > 1 and D(t 2 " ) = 1 otherwise. 
Our particular V has character 

ch(V) = £ t 2(i + ; + *)-2rf- 3) 

where 

B = {l,...,c}x {l,...,b} x 
The map D therefore yields 

, / / n / m n rr (/+/ + £ — d— l) 2 

_ C(c,b,a;d- l) 2 
C(c,b,a;d — 2) 2 ' 

where B' C B is the set of triples (i,j,k) such that 

i+7+it>d+3. 

Therefore 

C(c,b,d;d — 1) 



det(a(X)|_ 1 )=± 



C(c,b,d;d — 2) 



in thegiven basis. Combining this equation with equations I 
and (ra), we obtain 



N(a,b,c) = 



C(c,b,d;d — 1) 



C(c,b,d;d — 2)YlC(a,b;c — 1) 



Here, the symbol "fj" denotes repetition of the succeeding 
factor with cyclic permutation of a, £>, and c. 



We apply the inclusion-exclusion formula for each cube 
product and the inductive formula for the resulting simplex 
products: 



N(a,b,c) = 



T{3,d-l)I\T(3,a-2)l\T(2,a-l) 2 
T(3,d-2)~[\T(3,a-l)l\T(2,c-l) 
T(2,d-l)l\T(2,a-l) 
UT(2,c-l) • 



This is Propp's expression for N(a,b,c). 

Remark. Observe that a(X)\_ 1 is a diagonalizable, integer 
matrix with integer eigenvalues. The author considered a di- 
rect analysis of this matrix without using any representation 
theory other than the knowledge that the matrix can be de- 
composed using rational linear algebra. Unfortunately, no ob- 
vious pattern for the eigenvectors appeared in small examples. 
One way to express the coefficients of the eigenvectors with 
the aid of representation theory is to use Racah 3 j-symbols, 
which are complicated and have no obvious derivation using 
basic linear algebra. 

3. QUANTUM REPRESENTATION THEORY OF sl(2, C) 

The ordinary representation theory of sl(2,C) admits a q- 
analogue, the quantum representation theory, which can be 
used to ^-enumerate plane partitions. A proper algebraic treat- 
ment of the quantum representation theory would present it 
as the representation theory of a Hopf algebra f/ 9 (sl(2,C)) 
[S]. Briefly, an abstract representation category, or the rep- 
resentation theory of a Hopf algebra, is a category of vec- 
tor spaces on which some associative algebra or some matri- 
ces act, together with suitable definitions of trivial representa- 
tions, tensor-product representation, and dual representations. 
The most apparent difference between such a category and the 
ordinary representation theory of a Lie algebra is that A®B 
and B®A need not be isomorphic, and if they are isomorphic, 
they need not be canonically isomorphic. For brevity, we will 
avoid the definition of I/ ? (sl(2,C)) and give a direct, compu- 
tational definition of quantum representations of sl(2,C). 

Let He a complex number, and define q — exp(/i). More 
generally, for any number or matrix A, define q A = exp(M). 
Define the bracket of a number or matrix by 



[A} = 



g l/2_ 9 -l/2- 



The Laurent polynomial [n] is called a quantum integer and 
is related to the ^-integer (n) q by 



[n] = q 



A quantum representation of si (2, C) is defined to be a vec- 
tor space V and a function a : {H,X,Y} — * End(V) such that 

[a(H),a{X)]=2a(X), 
[a(H),a(Y)] = -2a(Y), 
[a(X),a(Y)] = [a(H)}. 



8 



The tensor product of two representation maps a and j3 is 
defined by the equations 

(a <g> P)(X) = a{X) ® qP W/ 4 + q - a WI 4 <g> j3(X), 
(a ® J3 ) (y) = a (F) ® ^ (H) / 4 + q - a WI 4 <g> p (Y ) , 
(a®P){H) = a{H)®I+I®${H). 

These equations are analyzed by Drinfel'd [0], but the dili- 
gent reader can check directly that the tensor product matrices 
again form a valid quantum representation. As a first step, 
[{a® P)(H)] involves q( a ®P)( H )/ 2 , and 

Ja®p)(H)/2 _ a(H)/2®I+I®P(H)/2 
q q 

_ q a(H)/2tSI q lesP(H)/2 
= q a WI 2 ®q^l\ 

since the exponential of a sum of commuting matrices is the 
product of the exponentials. 

The classical representation V„ deforms to a quantum rep- 
resentation by the formulas 



On(H) = 



OCn(X) = 



1 



Vo 

/0 [1] 
[2 







•• 


■ °\ 


n-2 


•• 


■ 





n-4 •• 


• 





•• 


.' - J 






[3] 








a„{Y) = 









■•• 


[n] 






Vo 





■•• 


oj 






(n 













°\ 


[n] 



















[n-l] 



















[n-2] 











\0 







[1] o/ 



Finally, although in the quantum representation theory there is 
no canonical isomorphism between V n <g> V k and V k (g) V n , they 
are isomorphic. Indeed, the Clebsch-Gordan theorem gener- 
alizes to the quantum representation theory. 
As before, let 



V = V 1 



-l 



6-1 



-H 



and define VU to be the eigenspace of tt(H) with eigenvalue 
A in this representation. In the tensor product basis, the ma- 
trix for a(X)\_ l is again a weighted adjacency matrix for 
Z(a,b,c). In the quantum case, the weighting has Kasteleyn 
curvature q at every hexagonal face. Therefore 

N(a,b,c) q = {deta(X)\_ l )/m g , 

where m„ is the weight of the matching corresponding to the 
empty plane partition. 



The rest of the derivation is a ^-analogue (or quantization) 
of Section \l.2\ The quantity m q is given by 



m„ 



> q = ±q-" bc/2 C(a,b;c- l) q C(a,c;b - l) q C(b,c;d- l) q , 
and the Clebsch-Gordan decomposition of V yields 



det(a(X)|_ 



± -abc/2 C(c,b,a;d-l) q 
C(c,b\a;d — 2) q 



We can q-ify the derivation at the end of Section 2.2 to prove 
the g-analogue of equation (Q): 



N(a,b,c) q = 



T{2,d- l) q T{2,a-l) q T(2,b-l) q T(2,c-l\ 
T(2,c-l) q T(2,b~l) q T(2,a-l) q 



4. OTHER DETERMINANTAL SYMMETRY CLASSES 



The method of Section 2.2 naturally extend to the enumer- 



ation of three other symmetry classes of plane partitions. 

Theorem 4.1. A Kasteleyn-flat, weighted adjacency ma- 
trix for each of the graphs Z KT (a,a,2b), Zp{a,a,a), 
and Zp KT (2a,2a,2a) arises as a Kt (X)\_ l , CCp(X)\_p and 
tt p KX (X)\_\ in representations CC KT , dtp, and CC p KT o/sl(2,C), 
with bases formed from weight bases of irreducible represen- 
tations. 

Recall the weighting of Z(a, a, 2b) induced by a(X)\_ l act- 
ing on V = V a+2b _ 1 <giV a+2b _i<giV 2a _ l . The symmetry ict 
preserves this weighting, and is induced by the linear trans- 
formation 



L KX e End(y„ 



+2*— 1 



■ V, 



given by 



a+2b- 1 1 



)Vr 



2a- 1 



k ' 



1, 



Define V KT to be the eigenspace of L KT with eigenvalue 
and choose vectors of the form e i ® e • ®e k — e ; - (g 

commutes with the action of 



) e, with 



i < j as a basis. Since L Kt 
sl(2,C), V Kt is also a representation. 

Recall that Z KZ (a,a,2b) is constructed by deleting vertices 
fixed by KT and edges incident to those vertices, and then tak- 
ing one of the two isomorphic components of the complement. 
By inspection of the basis for V KZ , the restriction a KT (X)\_ l 
of a(X)\_ l to V KX is a weighted, bipartite adjacency matrix 
for Z KT (a,a,2b). As in the no-symmetry case, opposite edges 
of each hexagonal face have the same weight, and therefore 
it is Kasteleyn-flat. Indeed, the weighting of Z KT (a,a,2b) is 
induced from that of Z(a,a,2b). It follows that 



N K% {a,a,2b)=tel(a KX {X)\_ l )/m % 



(8) 



where m KT is the value of a single term in the determinant (i.e., 
the weight of any matching in Z KT (a,a,2b)). 

To compute m KZ , observe that a fCT-invariant matching of 
Z(a,a,2b) has two edges of the same weight for each edge 
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in the corresponding matching of Z KX (a,a,2b), plus all of 
the deleted edges. The weights of the deleted edges are 
1,3,5, ... ,2a — 1. Therefore 



m KT = y/\m\/(2a-l)\\, 



(9) 



where «!! = («— l)(n — 3) . . . is the odd factorial. The de- 
terminant of (X (X ) | _ ] is best computed using the D map 
of Section 2.2. The representation V has two basis vectors 



e i ® e j <g> e k for each basis vector of V KT , plus basis vectors of 
the form e ( . ® e ; ® e^. Therefore 

i / a+2b 2a \ 

ch(V„) = 5 f ch(V)- g g^+2*-4«-4 6 -3j 



D(Pf^ - FT ft max(2; + fc-2 a -2fr-l,l) 2 
f=i iJmax(2/ + fe-l-2a-2Z7-2,l) 2 

_"+ 2i max(2;-2/j-l,l) 2 



max(2/ — 2a — 2b — 2, l) 2 
(2a + 2&-l)!! 2 
^-H^-l)! 2 ' 



Therefore 



D(ch(V KT )) 



l4 b - l {b-l)l 2 D(ch(V)) 



(2a + 2b-l)!! 2 
Combining with equations (§h and (pi), 



N KT {a,a,2b) 



lN{a,a,2b)2 b -\b-\)\{2a-\)\\ 



(2a + 2b-\)\ 



Define V p to be the eigenspace of L p with eigenvalue 1 , and 
choose as a basis for V p vectors of the form e t ® e ( - ® e ( - and 
vectors of the form 



when i, j, and k are not all equal. Since Z p (a,a,a) is the quo- 
tient graph Z(a,a,a)/p, the restriction a p (X)\_ l of a(X)|_j 
to Vp must be a weighted, bipartite adjacency matrix for 
Z p (a,a,a). Again, the weighting of Z p (a,a,a) induced from 
Z(a,a,a) is compatible with a p (X)\_ l and is Kasteleyn-fiat. 
If m p is the weight of any matching in Z p (a, a, a), then 



Let P(t) = ch(V) - 2ch(V Kt ) be the double summation; then Meanwhile, 



, , xx ,1, , ..3^max(6i — 3a— 1,1) 4 

Wp))=wch(v))n mm ; 6l -_ 3a _ 2 ; 1 i4 . 



because 



:h(y p )(0 = -(ch(y)( f )+2ch(y 2a _ 1 )( f 3 )) 



by counting basis vectors. The final result is 



/ " (3/-1) 2 

N p (a,a,a) = M(fl,fl,a)[] (3/_2)2 " 

This is again easily equivalent to the standard formulas jlO|]. 
In V = ^i 1 , we define 



This is equivalent to standard formulas for N KT (a,a,2b) 

Similarly, the symmetry p on Z(a,a,a) is induced by the 
linear transformation 



given by 



Lp (e f ® e, ® e t ) = e • ® e t ® e ; 



2a- 1 > 



together with the preferred basis of vectors e ( - A e • A e t with 
/ < 7 < fc. (Recall that the wedge product is the anti- 
symmetrized tensor product.) In this basis, a Kxp (X)\_ i is a 
Kasteleyn-fiat matrix for Z p KT (2a,2a,2a). A product formula 
for N p KX (2a, 2a, 2a) is the inevitable consequence, and in this 
final case, we omit the details. 
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